AD-H14 7  7 21  ELECTRONIC  STRUCTURE  CHEMICAL  BONDING  AND  ELECTRON  1/1 

CONDUCTIVITY  OF  THIN-F.  .  <U)  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  DEPT  OF  MATERIALS  SCIENC.  . 

UNCLASSIFIED  K  H  JOHNSON  ET  AL.  10  SEP  84  TR-8  F/G  20/12  NL 


A 147  721 


*CCU*STV  CLASSIFICATION  OF  T*IS  »AO(  r*hmt  Df  Ktttf  d) 


REPORT  DOCUMENTATION  PAGE 


4.  title  (m*  Mini*) 


Electronic  Structure.  Chemical  Bonding,  and 
Electron  Conductivity  of  Thln-Fllm  Transition- 
Metal  Slllcldes 


AUTHOMf.; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CAT  A  (.00  NUMBER 


i.  type  op  report  t  period  covcreo 

Interim 


4.  PERPORMINO  0RG.  REPORT  NUMBER 


OR  GRANT  NUMBERS 


K.  H.  Johnson  and  F.  A.  Leon 


Center  for  Materials  Science  and  Engineering, 
M.I.T. ,  Cambridge,  Massachusetts  02139 


II.  CONTROLLING  OPPICE  NAME  ANO  ADORES* 

Office  of  Naval  Research 
Department  of  the  Navy 


N00Q14-81-K-0499 


10.  RRQQRAM  (LtMCNT.  RROJfCT.  TASK 
A  RCA  A  WORK  UNIT  MUMlCRS 

Task  No.  Nr  056-757 


IS.  REPORT  OATB 

September  10,  1984 


IS.  NUMBER  fJ^PAGCS 


17.  OISTRIEUTION 


j^jgf  ,  Jfffimf  it  Miimmrnt  ham  Crnttrellku  0!H—>  IS.  SECURITY  CLASS.  (el  RW  J«P»rO 

Unclassified 


release; 

ted. 


STATEMENT  (mt  Ul«  ebec TMI  aniN  M  »!•«*  JO.  II  *«•"«  *«■  *«R«0 


is.  supplementary  notes 


Preuared  for  Publication  in  Solid  State  Communications 


I  (CeniNwe  mt  nntm  aide  It  iweeeaerf  I4mttlfk  kf  Hock  nunfrtr) 


electronic  structure 
chemical  bonding 
electron  conductivity 

thln-fllm  transition-metal  slllcldes  _ _ _ 


20.  ABSTRACT  (Cenltnue  an  nnm  bide  II  Identity  *P  *<•«*  number) 

^)The  local  electronic  densities  of  states  and  chemical  bonding  of  thln-fllm 
palladium  and  molybdenum  slllcldes  are  compared  on  the  basis  of  embedded  clus¬ 
ter  molecular-orbital  calculations.  Composite  Pd(d)-S1(p)ant1bond1ng/S1(p)- 
S1(p)bondlng  character  and  composite  Mo(d)rSi(p)nonbond1ng/Mo(d)-S1(d)bond1ng 
character  at  the  respective  Fermi  energies  are  responsible  for  the  different 
electrical  conductivities  of  these  slllclde  films  and  for  the  different 
Schottky  barriers  of  the  corresponding  slllclde/slllcon  systems.  — _ 


do 


1473 


COITION  OP  I  NOV  ts  IS  obsolete 
i  'H  01 02-LP-01 44401 


Unclassified 

SKCURlTr  CLAMIFICATION  OF  TwifRAOt  f"hmt  Dmm  Fnrered* 


a*'  '• 


>&/ 

OFFICE  OF  NAVAL  RESEARCH 
Contract  No.  N00014-81-K-0499 
Task  No.  Nr  056-757 


/£*  \\Licce55l° 

(\\  ♦  IKTIS  Gy 

\  •.  \  /  Ti'PT  n  'p  a  rj 


TECHNICAL  REPORT  NO.  8 


on  For 
A&I 

BTIC  TA3 
Unannounced 
Justif  icat  ion _ 


* 


By - - - 

_Dirtrihut ion/ 

Availability  Codes 
| Avail  and/or 
Spec ial 


Dlst 


l Ad 


ELECTRONIC  STRUCTURE,  CHEMICAL  BONDING,  AND  ELECTRON  CONDUCTIVITY 
OF  TH IN-FILM  TRANSITION-METAL  SILICIDES 

by 

K.  H.  Johnson  and  F.  A.  Leon 

Department  of  Materials  Science  and  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

September  10,  1984 

Reproduction  In  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


Approved  for  Public  Release:  Distribution  Unlimited 


ELECTRONIC  STRUCTURE,  CHEMICAL  BONDING,  AND  ELECTRON  CONDUCTIVITY 
OF  THIN-FILM  TRANSITION-METAL  SIL1CIDES* 
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Abstract 

The  local  electronic  densities  of  states  and  chemical  bonding  of 
thin-film  palladium  and  molybdenum  siltcfdes  are  compared  on  the 
basis  of  embedded  cluster  molecular-orbital  calculations.  Composite 
Pd(d)-Si (p)antibond1ng/S1 (p)-S1 (p)bonding  character  and  composite 
Mo(d)-S1 (p) nonbond 1ng/Mo(d) -SI (d) bonding  character  at  the  respective 
Fermi  energies  are  responsible  for  the  different  electrical 
conductivities  of  these  sllfcide  films  and  for  the  different  Schott ky 
barriers  of  the  correspond  1 ng  si  1 fcide/sl 1  icon  systems. 


Thin-film  transition-metal  slllcfdes  are  being  used  to  an 
increasing  extent  as  Interconnect  and  gate  materials  in 
very- 1 arge-sca 1 e- integrated  (VLSI)  circuits  [1-3].  Experimental 
[4-6]  and  theoretical  [7-9]  studies,  focussed  primarily  on  the 
neai — noble  Pd,  Ni ,  and  Pt/silicon  interfaces  and  silicides,  have 
established  the  importance  of  chemical  Intermixing  and  bond  formation 
in  determining  Schottky  barrier  height  and  electron  transport  across 
transition-metal  si  1 icide/si 1  icon  interfaces.  Such  issues  are  also 
expected  to  be  of  significance  in  refractory  transition-metal  (e.g. 
Mo,  W,  and  Nb)  si  1 icide/si 1  icon  interfaces,  which  are  of  great 
current  technological  interest  for  high-speed  VLSI  computer  circuits 
[10],  but  have  not  been  the  subjects  of  much  fundamental  experimental 
or  theoret i ca 1  i nvest i gat ion  [11]. 

In  this  communication,  we  present  quantitative  theoretical  models 
for  the  electronic  structures  and  chemical  bonding  of  thln-film 
palladium  and  molybdenum  silicides.  Because  of  the  importance  of 
local  chemical  bonding,  the  electronic  structures  have  been 
calculated  by  a  "real-space"  molecular-orbital  approach  rather  than  a 
"reciprocal-"  or  "k-space"  band-structure  method.  The  approach  used 
is  the  recently  developed  iterative  partitioned  scattered-wave  method 
[12],  which  permits  the  computation  of  molecular  orbitals  and  local 
densities  of  states  for  subclusters  of  an  extended  solid  or 
interface,  including  the  effects  of  "embedding"  the  subclusters  in 
the  extended  environment.  This  method  is  combined  with  the  Xalpha 
approximation  to  electron-electron  exchange  and  correlation  [13], 
which  has  traditionally  been  used  with  the  scattered-wave  method  in 
its  conventional  form  [14,15].  The  standard  Xalpha  scattered-wave 


molecular — orbital  method  has  already  been  applied  successfully  to 
elucidate  the  local  electronic  structures  of  hydrogenated  amorphous 
silicon  [16]  and  crystalline  silicon  containing  substitutional 
[17,18]  and  interstitial  [19]  transition-metal  impurities. 

In  the  iterative  partitioning  version  of  the  scattered-wave 
technique  [12],  the  secular  equation  for  the  molecular  orbitals  of 
large  silicide  aggregates,  such  as  the  Pd2SI  monolayer  thin  films 
shown  in  Fig.  1(a)  and  (b) ,  consisting  of  the  respective  encircled 
SlPd3  and  Pd^Si^  subclusters  and  the  surrounding  extended 
silicide  environment,  can  be  written  in  the  matrix  form: 


Here  T  represents  the  electron-wave  scattering  matrix  of  the 

O 

silicide  subcluster  and  Tb  represents  the  scattering  matrix  of  the 

surrounding  environment;  Gflb  and  are  matrices  which  represent 

the  "propagat i on"  of  electron  waves  between  the  silicide  subcluster 

and  its  environment;  and  A_  and  A.  are  the  electron  wavefunction 

a  d 

amplitudes  emanating  from  the  silicide  subcluster  and  environment, 
respectively.  To  make  the  solution  of  Eq.  (1)  practical,  we  first 
contract  the  matrix  of  (I)  using  a  simple  property  of  linear  algebra 
[20].  This  leads  to  the  matrix  equation: 


for  the  molecular  orbitals  of  the  sillclde  subcluster.  Includin' 


the  effects  of  coupling  to  the  sillclde  environment  through  the  terms 

G  .  T.  G. _ .  Solving  the  contracted  secular  equation  (2)  Is  much 

ab  b  ba 

more  computationally  efficient  than  solving  the  complete  secular 
equation  (1).  The  following  Iterative  procedure  Is  used  to  solve  Eq. 
(2):  first,  an  energy  eigenvalue  Eq  Is  found  from  the  subcluster 
submatrix  T  -1,  Including  only  the  molecular  potential  field  of 
the  surrounding  sillclde  environment.  Then  the  matrix  elements  of 
GabTbGba  are  ca1cu,at*d  using  EQ  as  the  trial  energy. 

Finally,  determinants  of  the  matrix 

tTa(E)]-'  -  Gab<E0)Tb<E0)Ob#(E0)  (3) 

are  calculated  at  various  values  of  the  energy  E  until  a  zero  Is 
found.  This  value  of  E  (call  It  Ej)  Is  assumed  to  be  an 
approximate  solution  of  (2).  If  Ej  differs  considerably  from  Eq, 
the  procedure  Is  repeated  from  the  second  step  using 


E  »  cEq  +  (1  -  c)Ej  (0  <  c  <  1)  (4) 

as  the  next  energy.  In  the  last  step  of  this  procedure,  the  matrix 

elements  of  G0b»  Tb,  and  are  unchanged,  so  they  must  be 

calculated  only  once  per  Iteration.  Therefore,  this  procedure  leads 
to  a  very  substantial  Improvement  In  computat I ona 1  efficiency  over 
the  direct  solution  of  Eq.  (1). 

The  resulting  local  electronic  densities  of  states  for  the 
encircled  subclusters  of  the  Pd9Sl  monolayers  of  Fig.  1(a)  and  (b). 


respectively.  Including  the  effects  of  "embedding”  the  subclusters  In 
the  surrounding  sllicide  environment,  are  shown  as  the  solid  and 
dashed  curves,  respectively  in  Fig.  2(a)  and  (b).  Consistent  with 
recent  band-structure  studies  of  epitaxial  Pd^SI  layers  on  a 
SI  (111)  surface  [9],  for  each  monolayer  there  Is  a  large  peak  in  the 
density  of  states  between  1  and  3  eV  below  the  Fermi  energy  Ep. 

The  peak  for  layer  (b)  of  Fig.  1  occurs  at  somewhat  lower  (more 
negative)  energy,  relative  to  the  Fermi  energy,  and  is  somewhat 
broader  than  the  peak  for  layer  (a)  of  Fig.  1.  This  is  due  to  the 
bonding  interaction  among  the  the  Pd(d)  orbitals  (1.  e.  "d-band" 
formation)  of  the  triangular  Pdg  clusters  in  layer  (b) .  The 
calculated  molecular  orbitals  corresponding  to  these  peaks  are  Pd(d) 
nonbonding  with  respect  to  the  SI (p)  orbitals.  The  secondary  peaks 
Immediately  below  and  above  the  nonbonding  peaks  correspond, 
respectively,  to  Pd(d)-S1 (p)  bonding  and  antibonding  molecular 
orbitals.  Indeed,  the  Fermi  level  coincides  with  Pd(d)-Si(p) 
anti bonding  states.  Within  layer  (a)  of  Fig  1,  these  anti bonding 

ft 

states  are  represented  by  the  SlPd^  subcluster  molecular-orbital 
wavefunctlon  contour  map  shown  In  Fig.  3(a).  Within  layer  (b)  of 
Fig.  1,  these  Pd(d)-SKp)  antibondtng  states  are  actually  bonding 
between  the  Si (p)  orbitals;  as  revealed  by  the  moleculai — orbital 
contour  map  shown  in  Fig  3(b).  In  other  words,  for  a  thin  PdgSI 
film  consisting  of  composite  layers  (a)  and  (b)  of  the  type  shown  In 
Fig.  1,  the  Fermi  level  Is  "pinned"  by  electronic  states  of  composite 
Pd(d)-S1 (plant ibondlng/SI (p)-Si (p)bondlng  character.  These  states 
lie  near  the  top  of  the  S1(p)-S1(p)  bonding  valence  band  of  silicon 
and  are  responsible  for  the  relatively  large  Schottky  barrier 
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observed  for  the  Pd2Si/S1  system  [11].  As  wl 1 1  be  demonstrated 
below,  this  composite  Pd(d)-Sf(p)  ant ibondfng/Sl (p)-Sl (p)  bonding 
character  at  the  Fermi  energy  also  largely  determines  the  electrical 
conductivity  of  a  Pd2Si  thin  film. 

The  calculated  local  density  of  states  for  an  SIMo^  subcluster 
of  an  MoSi2  monolayer  Is  shown  In  Fig.  4.  In  contrast  to  the  above 
results  for  Pd2Sl,  the  Fermi  energy  concldes  with  a  high  density  of 
Mo(d)  states  which  are  nonbonding  with  respect  to  the  Si (p)  orbitals. 
These  states  lie  well  above  the  Si(p)-SI(p)  bonding  states  that 
correspond  to  the  top  of  the  valence  band  of  silicon.  They  are 
responsible  for  pinning  the  Fermi  level  and  for  the  smaller  Schottky 
barrier  of  the  MoS^/Si  system  relative  to  that  of  the  Pd2S1/SI 
system  [11].  A  close  examination  of  the  moleculai — orbital  topology 
of  these  states,  as  exemplified  by  the  wavefunction  contour  maps  in 
Fig.  5,  reveals  that  the  Mold)  orbitals,  while  nonbonding  with 
respect  to  the  Si (p)  valence  orbitals,  are  actually  weakly  bonding 
with  respect  to  the  virtual  SI (d)  orbitals.  Although  one  usually 
Ignores  SI (d)  orbitals  as  being  relevant  to  solid-state  electronic 
structure-properties  relations.  In  this  situation  these  orbitals  act 
as  a  "pathway"  for  promot 1 ng  over  1 ap  and  meta 1 1 i c  bond i ng  among  the  d 
orbitals  of  the  Mo  atoms  which  are  only  second-nearest  neighbors  In 
the  sillcide  monolayer.  It  is  also  evident  from  these  results  (see 
below)  that  spatially  delocalized  Mo(d)-Si(d)  bonding  of  the  type 
shown  in  Fig.  5  Is  largely  responsible  for  the  electrical 
conductivity  of  an  MoS^  thin  film. 

The  computed  molecular  orbitals  can  be  used  as  a  basis  for 
calculating  the  electrical  conductivities  of  these  th in-film 
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suicides  via  Kubo  theory  [213*  Kubo's  formula  for  the  conductivity 
can  be  reduced  approximately  (for  "room  temperature")  to  the  form 
t22] 

2TI2ne2d2/h  (5) 

where  d  is  the  mo I  ecu 1 at — orbital  bond  distance  at  the  Fermi  energy 

and  n  is  the  bond  electron  density  at  the  Fermi  energy.  The 

Si(p)-Si(p)  bond  distance  in  Fig.  3(b)  is  3.8A,  whereas  the 

Mo(d)-Si (d)  bond  length  in  Fig.  5  is  2.6  A.  The  corresponding  values 

of  n.  obtained  from  the  computed  mol  ecu! at — orbital  components,  are 
22  22  -3 

3x10  and  2x10  cm  .respectively.  Substitution  of  these 

5 

values  of  d  and  n  In  Eq.  (5)  yields  conductivities  of  3.3x10  and 
4  - 1 

1.0x10  (ohm-cm)  for  the  PdgSi  and  MoSig  layers, 
respectively,  which  correspond  to  resistivities  of  30x1 06  and 
97x1 06  ohm-cm.  These  values  are  In  good  agreement  with  experiment 
[1].  Thus  the  difference  in  the  chemical  bonding  at  the  Fermi 
energies  of  PdgSt  and  MoSlg  thin  films  accounts  for  their 
respective  electrical  conductivities. 

I n  cone  1  us  1 on ,  it  has  been  shown  that  the  deta i 1 ed  compos i te 
mo 1  ecu 1  at — orbital  topology  of  the  electronic  states  around  the  Fermi 
energy  is  key  to  understanding  the  specific  electrical  properties  of 
thin-film  transition-metal  si  lie  ides.  Such  knowledge  may  possibly  be 
put  to  use  in  the  development  and  refinement  of  si  1 icide/si 1  icon  VLSI 


microstructures 
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Figure  Captions 


Fig.  I.  Monolayer  structures  of  PdgSI.  Solid  circles  represent 

Pd  atoms;  open  circles  represent  Si  atoms.  The  principal  subclusters 

|  defined  In  the  partitioned  scattered-wave  moleculai — orbital 

calculations  are  encircled  (dashed  circles). 

Fig.  2.  Subcluster  local  electronic  densities  of  states  for 

j  Pd2Si  monolayers  (a)  (solid  profile)  and  (b)  (dashed  profile)  of 

Fig.  1,  as  computed  by  the  partitioned  scattered-wave  method.  The 

principal  chemical  bonding,  nonbonding,  and  anti  bonding  characters  of 

^  the  component  molecular  orbitals  are  indicated. 

Fig.  3.  (a)  Contour  map  of  the  Pd(d)-St(p)  antibonding 

moleculai — orbital  wavefunctlon  at  the  Fermi  energy  for  the  SiPd3 

>  subcluster  of  Fig.  1(a)  monolayer;  (b)  contour  map  of  the  composite 

I 

Pd(d)-Si (p)anti bondlng/Si (p)-S1 (p)bondtng  moleculai — orbital 
wavefunctlon  at  the  Fermi  energy  for  the  Pd^SJ^  subcluster  of 
Fig.  1(b)  monolayer.  The  solid  and  dashed  contours  represent 

I 

positive  and  negative  values  of  the  wavefunctlon,  respectively. 

Fig.  4.  Subcluster  local  electronic  density  of  states  for  an 

MoSI-  monolayer,  as  computed  by  the  partitioned  scattered-wave 
I  c 

method.  The  principal  bonding,  nonbonding,  and  ant (bonding 
characters  of  the  component  molecular  orbitals  are  indicated. 

Fig.  5.  Contour  maps  of  the  SiMo-  subcluster  Mo(d)-S1(d)  bonding 

i 

moleculai — orbital  wavefunctlon  at  the  Fermi  energy  of  an  MoSI^ 

monolayer,  plotted  in  (a)  the  monolayer  plane  and  in  (b)  a  plane 

perpendicular  to  the  monolayer  and  containing  two  Mo  atoms.  Solid 

I 

and  dashed  contours  represent  positive  and  negative  values  of  the 
wavefunctlon,  respectively. 
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